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I Abstract. After a systematic introduction of some formulae for the energy radiated 

by localized electric charges and currents distributions, one considers the multipole 
, radiation and the reduction of the multipole tensors to the symmetric traceless ones. 

(D 
Q 

1. Introduction 

^ \ 

, Ph ! In the calculation of the energy radiated at large distances by a localised electric charged 
system it is not necessary to know the exact expressions of the electromagnetic fields E 
^ ! and B or of the potentials A and $. One may avoid the exact calculation, sometimes 
relatively complicate, in a simple way based on a formula for the power radiated by a 
charged system described by the charge p and current j densities with supports included 
in a finite domain V [Tj : 



dP, . ' " 



^ X ^^rad(^,^) 



(1) 



C/3 diZ /igC 

o ! 

. Here the origin O of the coordinates is chosen in the domain T), u = r/r, dP/dfi is 
^ ! related to the flow of the energy detected in the observation point r at large distance r 
Q^! compared with the dimensions of the given charged system. The vector Aj-^d is obtained 
from the retarded potential 

A(r.«)^£/i.V.«-f)dV, 

V 

with R = r — r', hj retaining only the dominant terms at large distances. A first 
approximation for this vector is obtained by retaining only the dominant term 1/r from 
the series expansion of l/R, 

= - + r' ■ V- +... = -- r' ■ V- + ... = - + + ... = - + 0(l/r2), 
K r \ RJ r'=o r r ^ j.^ ^ 

with a corresponding definition 

Aadir,t) = !^- [jir',t--)d'x'. (2) 
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In these relations and in the following ones we denote by 0(x") a series of powers of x 
beginning with x". 

Supposing r » X, where A is an arbitrary wave length from the radiation spectrum, 
such that the observation point is in the wave region, and retaining in equation only 
the terms having nonzero limits for r — > cxd, one obtains an approximate expression 
related to the energy flow observed in the point r and moment t. Rigorously, this is 
that part of the energy flowing in the neighbourhood of the observation point which 
contributes to the radiated energy. In the following we assume to work in this wave 
region. 



2. The radiation field 

In [ij the equation ((1} is justified using the supposed plane wave behaviour of the 
radiated field but also in P, in a footnote of the page 229, a rigorous proof is suggested 
for this. Indeed, this may be done by considering consistently only the terms from E 
and B contributing to the radiation j2]. 

Denoting by t' = t — R/c the retarded time, and by 

[p]=p(r',0, [j]=jir',t') 

the retarded charge and current densities, we have 



B{r,t) = V X A{r,t) = x 

4:71 



1 



r 

L V 



m'x' 



+ 0(l/r2) 



Air 



/io 1 
Att r 



[Vi)x /[j]dV-i/|[j]xVt'dV 

V V 

d 



+ 0{l/r' 



J{Vt')x-[j]d'x' + 0{l/r'). 

V 

Considering the series expansion of R, 

R = r + r' ■ {V'R)r'=o + ...=r-r'-Vr + ... = r 

we have 



r ■ r 



+ 0(l/r) 



t' = t- — = t-- + -iy-r' + 0(l/r), 

C CO 



and we may write 

Vt': 



1 r-r'l Ir' 
+ V- + + 

c c r c r 



1 



The vector B{r,t) may be written as 
J 

V 



dt 



+ 0(l/r) 



X 1/ +0(l/r2 



(3) 



(4) 
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so that the part of B contributing to the radiation is, in a first evaluation, 

1 / -^A^ad \ 



B 



rad 



X ly 



C\ dt 

The electric field E = — V$ — dA/dt with the retarded scalar potential 



(5) 



$(r,t) 



1 fp{r',t-R/cl^,^, 



R 



is given by 



E{r,t) 



. V 



£i/||,|dV-.0(l/.^) 



1 1 f9[p] ,3 / ^0 1 fd[j],3 ' , niM 2n 
' -ua x — : / d X +U{l/r ) 



Aneocr J dt 



V 



4:71 r J dt 



V 



where the equation Q is considered. Writing the continuity equation in the point r' at 
the retarded time t — R/c, 



|p(r'',t-^) + [V'-i(r',r)U. 



R/c — 0' 



and the relations 



V'j{r',t- 



R. 



[V'j(r', r)U,.n/c + " h{r\ t - ^) + 0(l/r) 



we have 



E{r,t) 



1 1 



c dV 
1 1 



' -a X + O^l/r ). 



An r J dt 



V 



The first integral in the right hand side of the last equation is zero because j = on 
the surface of V. Because {u ■ d[j]/dt)v — d[j]/dt = u x (v> x d[j]/dt), we obtain 



Eir,t) 



/io 1 

4tt r 



V X \ V X 



V 



d\j\ 
dt 



dV + 0(l/r'). 



As in the case of the magnetic field, we may write, in a first evaluation, the part of E 
contributing to the radiation: 

d_ 
di' 



(6) 



The equations verified by the fields Ej.^^^ and -Brad' 

~ ~ _ 1 _ _2 1 ~2 

-E'rad = C-Brad X i^, -Brad = X -Erad, ^O-^rad = — -^rad) 
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indicate the plane wave-like local structure of the radiation field corresponding to the 
radial direction u of propagation (these relations are valid also for the radiated fields 
Erad and Brad obtained from E and B by retaining only the terms with 1/r). 

Supposing that all the field variables are real functions, the Poynting vector of the 
radiated field ]^ ]^ 

-^rad = — -^rad X ^rad = ^O^lad = ~-^rad 
or ^ 

Hoc \ dt J 

Let the sphere of the radius r with the center in O. The radiated energy ^Al^i-ad passing 
through the surface element Aa, centered on the point r and corresponding to the solid 
angle AQ, in the time interval (t, t + 6t) is defined by 

6AW,^i = u-S,^ir'An6t. (7) 

From the last equation one sees that the angular distribution of the radiation power is 
given by the equation (jT]). 



3. The radiation of the point electric charge 

Usually, one derives the angular distribution of the power radiated by a point electric 
charge q using the results for the fields E and B obtained from the Lienard-Wiechert 
potentials and retaining from the corresponding expressions only the terms contributing 
to the radiation Here we illustrate the simplicity of the calculation using in this 

case the formula (jT}. 

For the sake of completeness we remind a concise introduction of the Lienard- 
Wiechert potentials Considering that the motion of the point charge q is given by 
the law t — y ^{t), the charge and current densities are represented as 

p{r,t) = q6[r-^{t)], j{r,t) = qv{t)6[r - ^{t)] 

where S is the Dirac function. The retarded vector potential A is represented by the 
integral 

^ ' ' A-K J \r-r'\ Att J J r - r'K 



hoo 



4:71 J R{t') 



Rit') 

t' -t + 



where R{t) = r — $,{t) gives the position of the observation point with respect to the 
particle at the moment t. In the last equation we may use the relation 
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supposing that the equation /(x) = has n roots Xi, . . . , x„ with /'(xj) = (d//dx) 
Here 

R{f) 



fit') =t'-t + 



(8) 



and 



fit') = 1 - m . nit') (9) 

where [3 = v/c and n(t) = Because f < c we have fit') > so that the 

equation fit') = has only one root r: 



1 



t-r--i?(r) = 0. 



(10) 



Denoting ?7(t) = 1 — /3(t) • nit), = riit)Rit), the vector potential A is given by the 
well-known Lienard-Wiechert expression 



An \sJt c 



0, 



and, obviously now, 



$(r,t) 



AttEq sir) 



Now let the potential Aj-^j corresponding to the field of the particle 



Aadi^, t) 



An r 



jir',t- R/c)d^x 



An r 



dt'vit')6 



t' -t + 



Rif) 



such that 



(12) 



An riiT)r 

with r defined by the equation (fTUj) . The time derivative of Aj.a,d is obtained by a simple 



calculation: 



9t 47rr 



a V dr] 

rj 7]'^ dt 



dr 
dt 



where a is the particle acceleration. To obtain dr/dt we consider the derivative of the 
equation (fTIH) 

dr dr (v ■ R\ 

di^ ^dt\ cR ) 



so that we get 



dr 
dt 



l-n-(3j^ rio' 



(13) 



where tiq = rjir). So 

dt Anr 



a V ( a ■ R v ■ n 



7] Tj^ 



cR cR 



1 

^0 
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and we may write 

^ = [r^a + (a . n)f3l + 0(l/r^). (14) 

at Arc rrjQ 

We have to calculate the expression v x dAj-g^f^/dt in order to introduce it in the equation 
(0). The observation point being given by its position vector r, it is clear that in the 
approximation considered here the same result is obtained for any chosen origin O inside 
the domain V. The unit vector i/ may be replaced by the unit vector n of the direction 
particle-observation point without changing the result for the angular distribution of 
the radiated power. This happens because 

n = ^ = iy-i^-V)- + ... = u-^ + ^i-^ + ... = u + 0(l/r) 
R r J. J. 



and 



h> X — - — = n X — V u[l r 

dt dt ^ ' 



n X 



By a straightforward calculation one obtains 

r]a^ + 2r]{n ■ a){f3 ■ a) - {1 - l3^){n ■ +0(l/r^)(15) 



Y dt J \ 4:71 J r'^riQ 

and with the equation ((T)) one gets the equation (73,9) given in pQ: 

g2 1 r 



dP 



r]a^ + 27]{n ■ a)(/3 ■ a) - (1 - /3^)(n • ay dfi (16) 



IGn'^eoc^ rj^ 

where all the particle parameters are considered at the retarded moment. Inserting the 
equation (fTHjl in the equation 0, and casting the terms of 0(l/r), we obtain that part 
of the energy passing through the surface element Act, in the time interval 6t, which 
contributes to the radiation. This is just the radiated part of the energy flowing in 
the solid angle AQ between the moments r and t. But this energy was emitted by the 
particle in the time interval St corresponding to the observation interval St. Therefore, 
writing 

dt 

SAW,^^ = 5rad ■ nr^An—Sr, 

the factor multiplying St in the right hand side of the last equation represents in fact the 
radiated part of the energy emitted by the particle in the unit of time (at an arbitrary 
time) in the direction n and in the solid angle AQ. Because dt/dr = r]{T) we may 
write the final result for the angular distribution of the radiation power emitted by the 
particle as an expression differing from the right hand side of the equation (fT^ by the 
factor l/rj^ instead of l/rj^ 



4. Multipolar expansion of Arad 

The series expansion of the integrand from the equation Q, by retaining only the 
1/r terms contributing to the radiation, leads finally to the multipole expansion of the 
radiation field. 
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Let us the Taylor series expansion of a function f{R), 



fiR) = E • ■■<A...^J(r) = E ^^^'1|V"/(r) (17) 



n=0 "-• n=0 



where 



■'il...«n 



doc 'i dx li 



and a" is the n-fold tensorial product (a0 . . . <2)ji...j„ = . . . ai„. Denoting by T*^"^ 
an nth order tensor, A^^^HB*^™) is an \n — m\th order tensor with the components 

By introducing the expansion (fTTj) into the equation we obtain 

n=0 



^rad{r,t) = £^eiE^^^n...*./<---<JiK^- ^)dV 



^°-E^^«S"^ (18) 



47rr ^0 n 



where 



(n) 



/<•••< J.K^-^)dV (19) 

and Si are the orthogonal unit vectors along the axes. 

In the following we use a generalisation to the dynamic case of a procedure given 
in [SI in the magnetostatic case. Let the identity 

V[xij{r,t)]=ji{r,t) + x,Vj{r,t). 
Considering the continuity equation + dp/dt = we may write 

Mr, t) = V[xd{r, t)] + x~pir, t) (20) 
and using this last equation in equation (fT^ .we get 



V V 

= ^:j(^',io)- V'«...OdV + 9,,...,„|<...<x:-p(r',to)d 



V V 

denoting to = t — r/c and considering a nul surface term because jt = on dV. Because 
of the symmetry of the derivative tensor and introducing the nth order electric multipole 
tensor 

p(«)(t) = j r"p(r,t)d^a;, (21) 

V 
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An) 



X. 



V 



V'^||-p("+i)(to) 



nd, 



j3 / 



V 



+ 

that is 



n + 



V«||^p("+i)(t ) 



(22) 



V 

By introducing the "vectorial product" T^"^ x a as the nth order tensor with components 



- ll...ln-li'^j 



and observing that, particularly, 



we may use in the equation (j22j) the definition of the nth order magnetic multipolar 
momentum given in j3] 



M(")(t) = X j{r,t)d^x. 

V 

So, the equation (j^^ may be written as 

= - V X (v"-^||M(")(to)) + ;7^^"II^P^"^'H^o) 

Going back to the expansion (fT^ . 



(23) 



(24) 



Anr 



VxY, V"-^||M(")(t - -) 

n=l 



_PU_ ^nn 9 p(n+l)/, _ ^^ 

47rr;^o(n + l)! "St ^ 



and, finally. 



47rr 



oo / 1 \n—l 



n=l 



47rr 1 n! 



i|p'"'(*-^) 



(25) 



Now, we extract from the equation (j25j) the terms contributing to the radiation. Because 

r (-1)" d" r 1 

--) = • • • - -) + O(-), 



Expressing the power radiated by electric charged systems 



we have in the first sum of the equation (j25|) : 



V X 



V"-i||M(")(t- -) 
c 



^i^'ijk(^j(^i\...in—i ^ii...in_ik (t 



-1)" 



d" r 1 



\rt-l 



n~l I 



\rt-l 



-U 



.n-1 1 



dt" ^ c 
d"M(")(t -r/c) 



X + 0(- 



dt" 



X u 



obtaining, finally, the formula 



Arad{r,t) = ^-y.^ 



d" 

"dt" 



X u 



47r r nlc'^ "dt" 



(26) 



giving the explicit contribution of each multipole to the radiation field. 

It is possible to express A^ad by the reduced multipolar tensors (total symmetric 
and traceless tensors) by applying the procedure given in [HI E! • This will be done in 
the following section. 



5. Expressing the radiation by reduced multipolar tensors 



For the following calculations it is suitable to write the expansion in the more 
explicit form 

1 



47rr _ 1 



n>l 



n\c" 



d" 



J k 



n>l 



i,-i||^pW(t ) 

" dt" ^ ^ 



d" 



n>l 



Vi- 



nyl 



n\c" 



_y. p. . 



(27) 



In ini Ej was done a general procedure for the reduction of multipole tensors represented 
by Cartesian components to symmetric traceless ones in the static and dynamic cases. 
The transformations implied by this reduction are defined such that the electromagnetic 
potentials A and $ are modified only by gauge transformations implying a specific 
feature of the dynanmic case: the redefinitions of the multipole tensors in the lower 
k < n orders induced by the reduction of tensors in a given order n. In the present 
paper this procedure is applied to the radiation field and, obviously, only the vector 
potential is to be considered. 

The reduction of multipole tensors begining with a given order n is achieved by the 
following steps. 

1. The reduction of the magnetic nth-order tensor M^^\ given by the equation 
to a symmetric tensor M(sym)- Since the magnetic tensor M*^"-* is symmetric only 
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in the first n — 1 indices, the reduction to a symmetric one may be performed by the 
transformation [H] 

1 

n 



^il...in ^ M {sym)ii .. .in 



Hi. ..I 



M 



11. ..I 



M, 



+ Mj^j2...i„_iii + • • • + Mj-^ 

-j^ n— 1 

A=l 

n-l 

A=l 

n—1 

- y M'^ ... - m(') . . 



A = l 

n-l 



ii...i„ 



n 



(28) 



A=l 



where we use the notations 



^il...in-l ^in-lpS^il...in-2PSJ 



fh...i„ fii...ix-ii\+i-in- (29) 

If M*^"^-* is given by the original definition ((221); ^ ~ Ith-order tensor N*^""^-* is given 
by 



N 



n 



n + 



n 



n + 



V 



(30) 



We write exphcitly the modification of the potential Aj-^fj induced by the substitution 



47rr 



A-KT Bi ""^ d" (A) 

A=l 



/io '""^"^ nine''' 



47rr 



'■rad 



n—1 jn 



nine 



A=l 



_ 71-1 jn 



47rr n — 1 

^rad + 



or 



where 



47rr 



^rad 



n!nc" 
47rr 



1 1 ^^„, V uy 



n — 1 
nine" 



An 

i."-i||^N("-i)(to) 



^rad + 



n — 1 



ij{r,t) 



liQ nine 
n — 1 



" dt" ^ ' 



4:Tcr n\nc^ 



d" 



" dr ^ ' 



+ ^ij{r,t)u (31) 



(32) 



and do not contribute to the fields ^^rad Bj-ad corresponding to a gauge 

transformation of the potential. 
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2. The extra-gauge alteration of the vector potential by the transformation (|31|) 
may be set off by the transformation of the electric multipolar tensor P^"^^); 

p{«-l) _^ p'("-l) = pin-l) _ d i^j(n-l) ^33^| 

c^n"^ dt 

such that the final transformation of the potential is the gauge transformation 

^rad ^ ^rad + (34) 

3. After the reduction of the magnetic tensor M^"-* to a symmetric one, we have to 

~(n) 

perform the reduction to a symmetric traceless tensor M . This reduction is achieved 
by the transformation [H] 

M(sym)ii...i„ Mil. ..in ~ M(sym)ii...i„ ~ ^ Siii2-^i3--in (35) 

Di^) 

where A*^"~^) is a symmetric tensor and the sum over D{i) is the sum over all 
permutations of the symbols ii,...,in which give distinct terms. Applequist |Hj has 
given an explicit formula for expressing the components of the symmetric traceless 
tensor M*^ ^ in terms of the traces of the tensor M {the detracer theorem, jH], equation 
(5.1)): 

Mil. ..in ~ l^(sym)ji...i„ 

K2W_iNm-l/2„ _ 1 _ OmV / ^ 

E l -L; y^n 1 z,nL).. ^ |v/|(n:m) 
/9„ _^\\\ 2^ "*1*2 ■ ■ ■ "i2m-li2m^^\sym)i2m+l-iu ^^^^ > 

m=l y^"- ^J-- D(i) 

where [n/2\ denotes the integer part oi njl and M|"y™^^,j^^^^ the components of the 
(n — 2m)th-order tensor obtained from M(sym) by the contractions of m pairs of symbols 
2, and obviously 

Er r _ \ \ " r r . «{n:'m) 

D{i) y /■ ■ P(i) 

(37) 

where the sum over P{i) is the sum over all the permutations of the symbols i. Using 
equation (IHUj) we may give explicitly the components of the tensor A*^""^-': 

[n/2] (_\Yi-\(1^ _ 1 _ 2rnV' , ^ 

m=l K'^"' i-)--'"' 

In terms of the tensor A the modification of Aj-^d induced by the substitutions 

and (jnni) is obtained by a straightforward calculation: 

Airr Attt Attt x a 

^rad ^rad r-r^iklJ^l^ii ■ ■ ■ ^i„-i ^iiia TT^^ia-in-ifc 

/io /io /io n!c" dt" 

47rr 47rr - 1 d"" . 

= ^rad + - ei——eikmPm, . . . Ui„_2-7-f^ii...in-2 

/io /io ^'c dt" 



(n- l)(n-2) d;^ 

2n!c^^ ^ikl^l'^h ■ ■ ■ '^in-3-^^h-in 

Anr Airr (n — l)(n — 2) 
Aa.d + + TTT^ ^iiki'i 



'dt" 
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So, the transformation of the potential may be written as 



Airr 



A 



rad 



A 



rad "I" 



Anr (n-l)(n-2) 
ihu + '-i -u X 



1/ 



n— 3 1 



(39) 



Ho fiQ /io 'znic"- ■ dt" 

4. It is a simple matter to see that the last extra-gauge term may be set off by the 
transformation 



|\/|{"-2) _^ |^'{n-2) _ |\y|{n-2) 



n - 2 d' 



(40) 



2c% dt2' 

5. This step consists in the reduction of the symmetric nth-order electric multipolar 
tensor P*-""-' to a symmetric and traceless one by a transformation of the type fl35j) : 



D{i) 



(41) 



where the symmetric tensor 11*^" is defined in terms of the traces of the tensor P^"-* 
by a relation similar to equation (jHHj) . The resulting transformation of Aj-^d is 



Anr Anr Anr 

Ho Ho Ho nlc"- ^ 



V', 



^'^ D(i) 



(42) 



In this equation we have (n — 1) terms with 5ik-, A; = 1, . . . , n — 1 and (n — l)(n — 2)/2 
terms with , j, = 1, 2, . . . , n — 1 so that 

47rr 47rr , Airr , n — 1 d" 

"^rad 



1^0 



Ho Ho 
(n-l)(n-2) 



2n!c 



,n-l 



-Ui, . . . Z/, 



n!c 

d^ 



that is 



where 



47rr 

fJ'O 



^rad 



Anr , Airr , (n — l)(n — 2) 



47rr (n — 1) 
Ho n!c"~^ 



1/ 



n-2 



|iLn("-2) 

'dt" 



2n!c' 



.n—l 



n— 3 I 



'dt' 



_j-|-(n-2) 



(43) 



(44) 



6. The alteration of the potential represented by the last term in equation ()43p is 
set off by the transformation 



p(n-2) _^ p(n-2) ^ ^ ~ n^"~^\ 



(45) 



2nc2 dt2" 

which preserves the symmetry proprties of P^"'"^). 

By this last transformation (|45p . the reduction of the multipolar tensors in the 
given nth-order is achieved. Now, to carry out this procedure to the n — 1th order, we 
must realize that in this order some tensors was been already modified in order to set 
off the alterations of the electromagnetic field by the reductions in the nth-order. So, 
the transformation ()33p alters the symmetry properties of the {n — l)th order electric 
multipole tensor because 

n — 1 d 



5p(n-l) 



c^n^ dt 



(46) 
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is symmetric only in the first n — 2 indices. To restore the full symmetry of the (n — l)th- 
order electric moment, we perform the reduction of N*-"'"^'' to a symmetric tensor by the 
transformation 

(n-2) 



Nji...i„_i — ^ Nji...j„_i 



N 



41. 



n 



1 

_ 1 ['^n...Jn-l'^li...l„_iJA 

^ A=l 



By introducing the tensor 7V^" '^^ with the components 
the transformation (|Trj)mav be written as 



(n-2) 



A=l 



and 



p: 



«l...«n-l 



P,: 
P 



n — 1 d 



n — \ d 



1 {"-2) J 
i^ti. ..«„_! ^ 2^2 ^^A^n-K? 1 .■'^^l...^„- 



n■■-n-l ^2^2 ^^--'l-'n-l ■ ^2^2 

If M^"-* is given by the original definition (j^^ . then we can write 

M.-.^^-. = / e'en • • • e.„-3 X d^'e 



25' 



n 



The alteration of the vector potential A by the transformation P9|l is given by 



A ^ A : U X 



nine" 



^-3^ ^^{n~2)^^^) 
U.6 



(47) 



(48) 



(49) 



(50) 



(51) 



This alteration of A is set off by the transformation of M given by equation ()4Up . 



|^'(n-2) ^ |\/|"("-2) ^ |\/|'("-2) _ 



n-2 d^ 



(52) 



n^(n — l)c^ dt^ 

By this transformation the symmetry properties of M are preserved. Particularly, 
by reducing the {n — 2)th-order multipolar tensors, in the case of M we have to 
achieve only the symmetrisation of the supplementary term from the equation (001) ■ 



6. Concluding remarks 

The equation ([Q) is a basic formula in the investigation of the radiation of electric 
charges distributions. In a course on Electrodynamics the radiation chapter is one of 
the mains goals of the introduction of Maxwell's equations. Sometimes it is necessary 
to give quickly some results illustrating the properties of the electromagnetic radiation. 
For the economy of this course it is a benefit to avoid some unnecessary intermediate 
results such as the expressions of the fields E and B. By exposing this problem in 
the sections 2 and 3, the present paper is not claimed as an original one but would be 
warranted at least by concerns of methodological completeness. 
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We summarize the results of the sections 4 and 5 by the following statements. 

(1) The reduction of the magnetic rath-order multipole tensor to a symmetric 
traceless one by the transformation (equations and (jSSI)) 

n—l 

M.,....„ ^ M,,...,„ = M,,...,„ ^^,^nM^U-^, " E ^ni^A^a-.i™ (53) 

^ A=l D(i) 

together the modifications of the electric {n — l)th-order and of the magnetic [n — 2)th- 
order multipole tensors, 

71 — 1 • 

where we use the super dot notation for the time derivatives and, 

n — 2 ■■ 

^ Mi,...,„_, + --^A,^.„,_^, (55) 

leads to a gauge transformation of the potential Aj-^d- 

We point out that if M*-"-* is given by the equation ()23p . the symmetric traceless 
tensor M may be identified with the tensor A^(") given by 

•^--W = („ + iK2i-i).. S/'-"" W-') " V)]..C,',.,.;d'x.(56) 

(2) The reduction of the electric nth-order multipole tensor to a symmetric traceless 

one , 

Pii...i„ ^ P = Pn...t„ - E ^^i^2ni3...i„, (57) 

together with the transformation of the electric {n — 2)th-order tensor, 

n — 2 ■■ 

Pii...i„_2 ^ Pii...i„_2 ~^ 2nC^ ^i\...in-2^ (58) 

leads also to a gauge transformation of the vector potential Aj-g^^j. 

We point out also that if P^"-* is given by the equation (j^H), the symmetric traceless 
tensor P*^ ■* may be identified with the tensor P^"-* given by jH] 

P..... = 7^^^ / P(r, t)r--^^ V"-d^x. (59) 

If we begin the reduction from a given order n, then the results of the reductions of P*^"-* 
and P^^-* are the tensors 7-*^"^ and A^^"-* given by the equations (|59p and (j56|) but for 
k < n the fcth-order reduced multipole tensors may differ from "P^^^ and Al^'^-' by terms 
iduced by the procedure of the reductions from the previous steps. These last terms 
give contributions to the potentials and fields expressed by toroidal moments and mean 
radii of various orders. 

We give here some simple examples of such reductions and we will see how naturally 
the toroidal moments appear as a result of such an approach. 

Let us the reduction of the magnetic and electric multipole tensors begins from the 
ixth and eth orders respectively (generally, considering the multipole's contributions of 
the same orders, ^ = e — 1). 
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For {fi,e) = (1,2), we have M^^^ ^ M^'^ = M«, P« ^ P^'^ = P«, P^^) ^ P^'^ = 
V^'^^ These transformations produce only a gauge transformation of A. 

For = (2,3), (3,4), (4,5), (5,6) the reductions are given in Appendix. For 

arbitrary /i and e, we think this is possible to find a general rule or, at least, to elaborate 
symbolic computer programs. 

In the case (/i, e) = (2, 3), 

p{i) ^ p(^) = p(i) - J_N^'^ + i-n(i). (60) 

4c2 6c2 

Here, N*-^^ and H^^^ are given by the equation ()A.10|) for Pg^j = 0, Ng^j = 0, Pqqppi = 
that is eliminating the contributions from the orders > 2 of the magnetic multipole 
tensors and from the orders > 3 for the electric ones. 

Taking into account the continuity equation verified by p and j, we obtain 

P^ = P. - (61) 



c2 



where 

[«-j){.-2«'j.ld'?, (62) 
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V 

is the toroid dipole tensor [10-12]. 

In the case (/x, e) = (3, 4) we have the changes 

p(4)^ p(3) p(3) ^ p(3)^ p(2) ^ p(2) ^ p(2) _ An^^^ + — n(2) 

9c2 4c2 

P« - An^'^ + Aii^^^ (63) 
4c^ 6c^ 

N., = ^(N., + N,,) (64) 

and Aj, Njj and Iljj are given by the equations ()A.3|) . ()A.4|) . ()A.7|) and ()A.8|) by 
eliminating the contributions from the orders > 3 and rig > 4. In this case one 
obtains the contribution of the toroidal quadrupol tensor T^'^\ [10-13]: 







p(4) _ 


p(^) 


p(l) ^ 


p(^) 


where 





42 



^ / [4(^ • JMk - 5ei^^Jk + ikji) + 2e{i ■ Mk] d'^ (65) 



V 



having, beside the equation (pT|) . 



Pik = Vik (66) 



c2 



and the dipolar magnetic moment modified by a mean-square current radius: 

1 1 



M. + ^-/e^(^x,-)d^ 

V 
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In the case (yU, e) = (4, 5) we obtain the following results of the reductions: 



M 



(4) 



(2) 



M 

p(3) 
S(2) 



42c2 
1 



/ e [e. X j), + X j) J d% 

V 



20c2 



(67) 



V 



p(5) ^ p(4) 



p(3) 
1 

60 



3 ~(3) 3 ~(3) ~ 
-N + : P,,fc 



V, 



ijk 



V 



16c2 10c2 " " 

D{i,j,k) D{i,j,k) D{i,j,k) 

1 



P 

A 



(1) 



2 ~(2) 1 ~(2) 

4c2 96c4 



P. 



-p. . T . . 

It] (.2 *J ' 

■p.-It.-Ia- 



7? 



These results show that one may obtain from the formula ()26|) the correct 
representation of the electromagnetic field by the reduced multipolar tensors but 
introducing these tensors up to a given order n, we obtain separate contributions from 
some electric toroidal moments and mean 2?2— power radii. This was pointed out firstly 
by Dubovik et al [ 10-12]. In the present paper we point out that considering the 
contributions to the electromagnetic field of some toroidal moments, one suppose the 
reduction of the multipole tensors up to a well defined maximal order n. 
We illustrate also this statement by calculating the total power radiated by a system of 
electric charges and currents. 

Let us the total radiation power obtained by integrating the equation (^: 



X, 



1 



J/ X A) rMn(i/) 



(6^ 



considering only the contributions of the magnetic and electric multipoles up to the /ith 

and eth orders respectively. Usig equation ()26p we may write 
2 



/ 47rr 
V /^o / 

= EE 

n=l m=l 

e £ 

+ EE 

n=l m=l 

M £ 



[u X A) 



I fi,e 



?T,!m!c"+'" 
1 

l™l^n+m-2 



(m) 
m 



n\m\c 



(m) 
m 



(69) 



n=l m,=l 
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where 

(n) _ d 

The calculation of the integrals in equation is reduced to the calculation of 
< i^h--- i^ir, n = 0, 1, . . . with 



< I/jj . . . Ui^^ >v— Cn ^ (^iiJa ■ ■ ■ '^«2n-li2n5 ~ (9 _L 1 "ill 

Let us the symmetric traceless tensors A*-"-* and B*-'"-' and the averaged contraction 



(70) 



jy'^IIA^")') II fjy'^'lIB^"') 



J^il ■ ■ ■ ^ik^jl ■ ■ ■ ^jkl^il-..ikik + \---inBji...jyjy^-^...j,n/^ ■ 



This is non zero only for the products of with p = 1, . . . k, q = 1, . . .k' and it is 

k\ 



easy to demonstrate the relation 

'(i/*^||A("))||(i.'='||B(™) 



/\(")||gM 



^k'k- 



(71) 



V (2A;+1)!! 

The terms of the last sum from the equation ()69|) give contributions to the total radiated 
power of the form 



■ ■Jm-iq 



but all the terms from the sum of 5— products representing the averaged products of 
z/'s contains either 6i^p or 6pj^, k = 1, . . . ,n — 1, I = 1, . . . ,m — 1 such that, because of 
Ei^pq and of the traceless character of A and B, the result is zero. Using these results in 
equations (jHHj) and we obtain 

(72) 



X, 



47reoC^ 



E 



n + l 



{nn\{2n + l)\\c 



2n 



'^',n+lll'V',n+l 



+ E 



n + 1 



\nn\{2n + l)\\c 



2n-2 



pin) p(n) ■ 



For comparison with results existing in literature [1,10-12] we write here the results in 
the following cases. The case (/x, e) = (1,2) is given in jj: 



2^1,2 



1 



-P2 P2 + — PV + MV 



1 ...(2) ,, ...(2)' 

V \\v 



20c2 



(73) 



this result being justified by the invariance of the radiation field to the transformation 

p(2) ^ p(2)^ 

In the case (/x, e) = (2, 3) we obtain 



2,3 
1 



2K(i)Mf^a) , 1 K(2)mK(2) 



ttPVIIPV + 



20c2 



+ 



^(3)mK(3) 



1 



7- 4 iiPVH ^Mo ||Mn H -,v,o 

945^4 ,4 11 ,4 ^ ,2 11 ,2 ^ 20c4 '3 



(2) ,,^(2)' 



M\'||M 
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-\P 



+ 20?-^ 

In the case (/i, e 



1 ■■■ 2 



945c4 
(3,4), 



3c2 

2 p(3)||^(3) 



20c2 



V 



(2) 



(2) 



2^3,4 
1 



3 2 11^-2 +20c2'^'3 11^^,3 + 945^4- ,4 



2 ~(3),,~{3) 

PVIIPV + 



18144c6 



>W||p(4) 



1 



3V^-^ 
1 



T 



20 r C2 dt4 



2 /.. 1 V 

+3?r+?d^^j 

' r c2 dt4 r 20 



945c4 

1 



_L I _p(4) 1 1 ±I_p(4) 

18144 c6dt5 "dt5 



M 



(2) 



(74) 



(75) 



Appendix A. Reduction of multipole tensors 

In tis Appendix we give in a diagram form the reductions of multipole tensors for the 
cases (/i,£) = (2,3), (3,4), (4,5), (5,6). 



(/i, e) = (2, 3) 



p(3) 



n 



(1) 



p{3) 
p(2) 

nl 



Eq.J33,n=2 



p(i) ^ p(i) _ _1^N 



(1) 



Eq.0nj."=3 



p{2) 
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(/i, e) = (3, 4) 



'V"(sym) 



A 



(1) 



N<4 



p(4) 



n 



(2) 



p(4) 



p(3) 



n 



(1) 



p{3) 



Eq.0SJ,n=3 



Ea.l^.n=?, 



p(2) ^ p(2) _ 



(2) 



EQ.i^.n=2 



p(l) ^ p(l) _ 1 |\|(^) 

4c2 



Ea.BSt.n=4 



Eq.0SJ,n=3 



p(i) _ A^U'" 



p(2) _ 


9c2 ^ 4c2 


n(2) 




p(2) _ 


9^"sym 1 4c- 


.n(2) 



Ea.<52l.n=8 



n 
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(/i, e) = (4, 5) 









MS.) 




i 



M(2) + ^A(2) 



'^'(sym) + 4c2^^ 



p(5) 



n 



(3) 



-p(5) 



p(4) 



n 



(2) 



Ea.i33l.n=4 



n=4 



Ea.l^.n=3 



Ea.W]i.n=3 



Ea.i33l.ra=2 



Ea.E5l.n=5 



p(3) , p(3) 3_M 



(3) 



M(2) ^ M(2) + ^A(2) 



p(2) ^ p(2) _ 



(2) 



(1) 



m(3) 



Ea.i45l.ra=4 





p(4) 








1 


p(3) 3_m(3) _3_fr{3) 

^ 16c2 '^sym ^ 10c2^^ 







Ea.i52l.ra=4 



Ea.B^.n=3 



p« - ^N^^^ p« - ^N^^^ + ^n^^) 
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'(1) 



n 



N 



(1) 



~(2) 



Ea.<^.n=3 



M«.. 



EQ.rOl.n=2 



p(l)^ 
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= (5,6) 



Eq.03,n=5 



Eo.Bni.n=5 



Eq.E3,n=4 



M 



(4) 

(sym) 



n=4 











|U|(3) I ^A(3 






A<(3) + 


, ^(3) 










'V"(sym) + 4c2^^ 



Eq.J23,n=3 



n=3 



Ea.i33l.n=2 



p(6) 



n 



(4) 



Ea.i45l.n=6 



p(4) 



(4) 



p(3) ^ p(3) _ 



M(2) + 1 A(2) 



p(2) 



M(i) + 1 Ad) 



(1) 



p{6) 



p(5) 



n 



(3) 



Ea.BSl.n=5 



7?(5) 



Expressing the power radiated by electric charged systems 



23 



p(4) _ 4 ^(4) + 1 n{4) 



p(4) _ ^m(^) , ^n(4) 



p(3) _ 3 , 3 n(3) 

^ 16c2 '^sym ^ 10c2^^ 



n 



(1) 



, ~(3) , ~(3) 



n 



+ 10c2^^ 100c2-^^ 



N 



A//(3) I _3_ a' 3_yW-(3) 

+ 10c2^^ 100c2-'^sym 



Ea.l5^.n=5 



Ea.lISl.n=4 



Ea.lE^.n=4 



Eq.E3,n=3 



Ea.lS^.n=3 



A 



Eq.OSJ 



p(2) _ ^N^'^ ^ p(2) - 2 + 1 n(2) 



9c2 



P« - 



(1) 



(1) 



N 



,'(2) 



IVI .. — > IVI -h g^a^i ^g^a^V 200c4^^ 



~(3) ~(3) 



A^(2) + _i^A(2) _ 



N 



Eq.OSJ 



p(l)^ 



~(2) 

AW (2) I J_ a' ^AT^^) 

•/^l ^ 4c2^^ 24c2-'^sym 
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o 1 ~(2) -, -'(2) 



-P 2) 2_|U , J_TT J V^-l 

In these diagrams are used the following notations: 

V V 



V 

D{i) ^ D(i) 



, 1 

^J 2° 



' qqppk 

D{i) 



V 



3 •• 

—M, 



D(i) ^ 



3^ ..xi^X3)]id% 

V 



nit) 



n -^p - ^ 

'-Hjk g qqijk 



+ 



D(i) 

Q X/ ^ij^ qqppk: -^ijk ~ ~'q J ^ 

D(j) J) 

l^^ijPqqpp - (j^ii^j + l^ggii) 



[Ip - -Up 

Y qq^J YQ ^^^^ 
1 rl" 1 ■• ■ 

Y^iqij ~ j^^ij^qqpp 



^ik--\j i% a X 3)k d'e, J^i^-\je{ix 3\ d% n,,- = n,,. - ^-5 



V 

V V 
V 

and, 



A - A - 1 



^ D{i) 

\sym)ijkl - E {^qqkl + ^qqlk) , 



^ijkl — ^ijkl — E [ 7^55^; " E ^^^^5^^^ 
D({) V 0(0 
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^ijk ^(syni)ijk r ^ ^ ^ij^qqky 

D{i) 

^ijk ^ijk T ^ ' ^ij^qqk) 
^ D(i) 



J^ijk A/'(sym)ijfc ^ ^ ^ ^ij-^qqk- 

Dii) 

For /i < 5 and £ < 6, the quantities N, A, ... are obtained from the above expressions by 
ehminating the contributions of the magnetic multipolar tensors of orders Um > and 
of the the electric ones for orders Ug > e. 

We point out that the reduction diagrams described in this appendix are valid also 
in the case of an arbitrary electromagnetic field as may be seen from [Zj. 
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